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Abstract. In this paper, we prove the off-shell equation satisfied by the transfer matrix associated with the 
XXZ spin-i chain on the segment with two generic integrable boundaries acting on the Bethe vector. The 
essential step is to prove that the expression of the action of a modified creation operator on the Bethe vector 
has an off-shell structure which results in an inhomogeneous term in the eigenvalues and Bethe equations of 
the corresponding transfer matrix. 
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1. Introduction 

In the last decade, the study of the spectral problems associated with quantum integrable models which 
are not invariant by U{ 1) symmetry has been an active research field. The prototype model is the XXZ spin- 
i chain on the segment with general boundary conditions, see II El for the Hamiltonian in our notations. 
For non-diagonal boundaries, although the model is known to be integrable, the Bethe ansatz techniques, 
and in particular the algebraic Bethe ansatz (ABA) [29:, [30], could not be directly applied due to the lack 
of a reference state for the transfer matrijfl. There were many attempts to overcome this problem: by 
introducing constraints between the left and the right boundary parameters 0 m ES EE by considering 
specific non-diagonal integrable boundaries^ and introducing a new structure for the Bethe vector [5] [21], 
or by developing alternative methods to the Bethe ansatz Him nu, including the separation of variables 
(SoV) [26] [16] method. More historic details can be found in [4 . 

Recently, an essential step to construct the solution of quantum integrable models without 17(1) symmetry 
by means of the Bethe ansatz was accomplished by the introduction of the so-called off-diagonal Bethe ansatz 
(ODBA) [UJ [XU]. The key point is the use, as an ansatz to characterize the spectrum of such models, of a 
modified T-Q Baxter equation which exhibits an additional inhomogeneous term. 

These results brought the questions: what are the Bethe vector associated with the new eigenvalues and, 
in addition, can the new inhomogeneous term be obtained solely from an ABA perspective? 

The first progress arose in [5], where the Bethe vector of the XXX spin-^ chain on the segment were 
conjectured from the modified T-Q Baxter equation [10]. In HE] the origin of the additional inhomogeneous 
term was clarified by means of a modified algebraic Bethe ansatz (MABA). In the MABA we drop the quest of 
the reference state and just focus on the construction of a highest weight vector for the underlying reflection 
algebra. In this framework, the additional term is understood as the action of a modified creation operator 
on the off-shell Bethe vectoi0. Such off-shell action was proved for the totally asymmetric simple exclusion 
process (TASEP) model (a degenerate case of the XXZ spin-^ chain on the segment with non-diagonal 
boundaries) in [T2J. This proof uses the so-called F-basis 22] to analyze the pole structure of the creation 

1 For models with U( 1) symmetry and q not root of unit, the reference state corresponds to the highest weight vector of the 
finite dimensional representation of the underlying algebra. 

2 In a way such that a reference state for the transfer matrix can still be identified. 

3 A Bethe vector is said to be off-shell if it depends on arbitrary parameters and on-shell if these parameters satisfy the 
Bethe equations. 
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operator. Let us also mention that a proof of the Bethe vector expression, solely for the on-shell case, was 
proposed from the ODBA for the XXX, in [33] . and XXZ, in [34], spin-^ chains on the segment. In this 
case, the canonical SoV-basis m was used. 

Here we prove the validity of MABA approach for the transfer matrix of the finite XXZ spin-i chain on 
the segment with two generic boundaries. The expression of the off-shell Bethe vector of this model were 
conjectured in [7]. Let us recall that the construction of the Bethe vector follows the steps: 

• Find the highest weight representation for the operators that describe the model through the general¬ 
ized quantum inverse scattering method [30] . This was done in [8] from a local gauge transformation. 
This is possible if the dynamical right reflection matrix is diagonal or upper-triangular. 

• Use the gauge transformation to map the transfer matrix to a modified lower triangular form [7] or 
to a modified diagonal form [5j. This step introduces the modified operator used to construct the 
Bethe vector (which is built from applying successively the modified operator on the highest weight 
vector). 

• Apply the transfer matrix on the Bethe vector. Use the commutation relations and the action of the 
modified operators on the highest weight representation to obtain the off-shell action of the transfer 
matrix. This action has a new term with an additional modified creation operator acting on the 
Bethe vector. 

• Show that the new term has a wanted-unwanted structure. From direct resolution of chains with 
small length (N = 1,2), we notice that this is possible if the number of creation operators equals the 
length of the chain. We observe that these wanted-unwanted terms give the inhomogeneous part in 
the eigenvalues and in the Bethe equations of the transfer matrix introduced in the context of the 
ODBA. This step was performed in [7]. 

In this paper, we prove the explicit form of the off-shell action of the modified creation operator on the 
off-shell Bethe vector using the canonical SoV-basis m- As an additional result, we make the explicit link 
between the MABA and the SoV [16) [26] and ODBA jlOl [331 S3] methods. 

The paper is organized as follows. In section [2] we recall the definition of the transfer matrix of XXZ 
spin-i chain on the segment in the generalized quantum inverse scattering framework as well as the gauge 
transformation that yields the dynamical operators. In section[3] we recall the highest weight representation 
of the dynamical operators and introduce the modified dynamical operators. In section [4] we implement the 
MABA and in section [5] we prove the action formula of the modified creation operator on the off-shell Bethe 
vector. The connection with the ODBA and SoV methods is given in section [6] We conclude in section [7] 


2. Generalized quantum inverse scattering and gauge tranformation 

The transfer matrix of the finite XXZ spin-i chain on the segment is constructed following the generalized 
quantum inverse scattering framework [30] . Here we recall only the essential features, more details in our 
notations can be found in [HE]- From the trigonometric I?—matrix, which acts on V a ® vE 


( 2 . 1 ) 


A,:,,?, (w) 


( b(qu) 0 0 0 \ 

0 b(u ) 1 0 

0 1 b(u) 0 

\ 0 0 0 b{qu) / 


b(u) 


u — u 

_ 1 1 

q-q 


^In this work, we only use two-dimensional complex vector spaces, i.e. V = C 2 . 
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and from the K —matrices HE which act on V, 

(2.2) K~(u) = 

(2.3) K+(u) = 


k ( u) t 2 c{u) \ , — / . _! 

.2 1) ’ k {u)=v-u + v + u 


T 2 c[u) k (u *) 
k + (qu) R 2 c(qu) 


' c(qu) k + (q 1 u 1 ) 


c(u ) = u 2 — u 
, k + (u) = e + u + e-U ^ 1 , 


where {e±,K, k} and {p±,t, f} are generic parameters, one constructs the transfer matrix, that acts on the 
Hilbert space R = 


(2.4) 

where 


t(u) = tr a (K+(u)K a {u)), 


(2.5) 

( 2 . 6 ) 


Ka{u ) = R a \{u/vi)... R aN (u/v N )K a ( u)R aN {uv N )... f? a i(m;i), 


£/ ( u) &(u) 

V(u) @{u)+ 5 ^fy^(w) 


is the double row monodromy matrix that acts on V a 0 R. This formalism allows one to reformulate the 
problem in the operator language. The transfer matrix is given by, 

(2.7) t(u) = (f>(u)k + (u)£ /( u ) + fc + (g -1 w -1 )i^(ii) + c(qu) ^ k 2 &{u) + R 2 ^(u)^ , 


with 

( 2 . 8 ) 


<t)(u) 


b{q 2 u 2 ) 

b(qu 2 ) 


An important point in the ABA framework is the existence of a highest weight vector which is annihilated 
by the operator ^( u ) and which is an eigenstate of the operators 32 /( u ) and 3>{u). This vector exists when 
we have an upper triangular K~ matrix. However, for generic K~ matrix, one needs to use the gauge 
transformation mi in order to define new double-row operators that are dynamical, i.e. they depend of an 
additional parameter said dynamical, and to be able to construct such highest weight vector. In addition, 
the gauge transformation allow us to eliminate the dynamical creation and/or annihilation operators in the 
transfer matrix expression. The building blocks of the gauge transformation are the covariant vectors, 


(2.9) 


I X(u,m)) = 


aq 


-m n , — 1 


I Y(u,m)) = 


/V 


and the contravariant vectors, 

(2.10) (X(u,m) | = —^—— ( —1, ), 

7m-l 


(Y(u , to) | 


qu 

7m+l 


1, 


—Pq m u 1 ) , 


where a and /3 are generic complex parameters, to is an integer and 7(11, to) = aq m u — /3q m u 1 with 
7(1, to) = 7 m . They satisfy the shifted scalar products, 

(2.11) (X(u,m)\X(u,m)) = (Y(u, m)\Y{u, to)) = 0, 

(2.12) (X(u,m+ 1)| Y(u,m — 1)) = ( Y(u,m — 1)|X(m,to+ 1)) = 1, 


as well as the shifted closure relation, 

(2.13) |Y(u, to — l))(X(u, to + 1)| + \X(u, to + 1))(F(m, to — 1)| = ^ ^ l)' 
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From these vectors we can introduce the dynamical operators, 

(2.14) < tf(u,m) = (X(u, m)\K a (u)\X(u~ 1 , in)), 

(2.15) &(u,m) = (Y(u,m)\K a (u)\Y(u-\m)), 

(2.16) £/(u,m) = (Y(u,m — 2)\K a (u)\X(u -1 , to)), 

(2.17) @(u,m) = ^ m+1 (X(u,m + 2)\K a (u)\Y(u~ 1 ,m)} - ^ £/(u,m). 

7m b{qu 2 )'j m 

such that the transfer matrix can be written in the following way, 

(2.18) t(u) = t d (u , to) + u~ 1 c{qu) (c, m £g{u, to) - ( m &(u, to)) , 


with 

(2.19) ( m = -(a9- ra - 1 +4)(a(i- ra - 1 +4), Cm = — (ftT" 1 + &) (ftT" 1 
and the new parametrization, 

(2.20) e_ = iKK(€/£+£/£), e+=*««(££ + 1/(CC)). 

The term t d (u,m), which is the diagonal part of the transfer matrix, is given by 



(2.21) td(u, to) = a(u, m)jz/(u, to) + d(u, m)0(u, to), 
where 

(2.22) a(u,m) = a(u) — 5 m (l){q^ 1 u~ 1 )u^ 1 c{qu), d{u,m ) = d(u) + u~ 1 c(qu)S m 
with 

(2.23) a(u) = u -1 </>(w)/c + (u), d(u) = u _1 fc + ((7 _1 « _1 ), 

(2.24) k + (u ) = ihn{^u + £ _1 w _1 )(C _1 m + Cw -1 ), 
and 

(2.25) <5 m = — (ag-™- 1 + (/?g m+1 +i^). 

7m+i \ k £' v 

Let us introduce finally a bookkeeping notation. For generic gauge parameters {«,/?} the dynamical 
operators are denoted by, 


{ 42 /(u, to), &(u, to), ^(u, to), ^(u, to)}. 

If we consider some fixed gauge parameters {ax,/dr} we denote the corresponding fixed dynamical operators 
with indices, namely, 


3. Representation theory and modified operators 

In order to consider the diagonalization of the transfer matrix in the framework of the generalized quantum 
inverse scattering method, we want first to introduce a suitable basis to construct the Bethe vector. In 
particular, since we want a highest weight representation, it can be built for the dynamical operators with 
the choice of gauge parameters 

a^a dr = iq m ° +N ^, /3 ^/3 dr = iq~ m °~ N ^, 

T fl Tfl 


(3.1) 
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which corresponds to bring the dynamical right K -matrix to a diagonal form, see (8j for details. As a 
consequence, it allows us to construct the highest weight vector, 

(3-2) l^m 0 )dr =®iL 1 \X(Vi,mo+i))dr- 

such that the actions of the dynamical operators at the point m = mo are given by, 


£&dr ( w , Too )I) dr — [u)A[u)\tt mo )dr, 
3>dr(u,mo)\n")dr =U(j)(q~ 1 u~ 1 )k~(q~ 1 u- 1 )h{q~ 


m o/ 

"(it, mo)\£l mo )dr = 0, 


L )I0 


dr 5 


N 


A(u) = W_b(qu/vi)b{quVi), 


(3.3) 

(3.4) 

(3.5) 
where 

(3.6) 
and 

(3.7) k~(u) = iTT(fj 1 u + + flu- 1 ), 

with the new parametrization 

(3.8) = *fr(/x//2 + jl/fj), v+ = ifT^fxp, + 1/ (/z/i)). 

Moreover the dynamical operator 3§dr is dynamically nilpotent of order TV + 1, 

(3.9) 3§dr(u,m 0 + 2N)£S dr (ui,m 0 + 2(TV - 1 ))... & d r(u N ,m 0 ) = 0, 

where u = {rt 1; ..., mat} is a set of TV generic parameters. 

Let us now consider the string of M operators acting on the highest weight vector, 

(3.10) @dr{uh,m 0 + 2(M - 1)) .. .&dr(Uj M ,mo)\ft™ 0 ) dr 

where M £ {0,, TV} and {uj 1 ,..., Uj M } is some subset of cardinality M of u. Since the vectors (13.101) are 
independent for generic u, from these vectors and for fixed M we generate a subspace Wj^ 0 of dimension 

^ ^ ^ Hilbert space. The highest weight representation induces the decomposition 


N 

(3.11) H=0w;° 

M—0 

of the Hilbert space. Indeed the total dimension of the sum on the r.h.s. is 

< 3 - 12 > E(m)= 2 " 

M=0 

Thus, in order to construct the Bethe vector of the model, we are looking for a linear combination of the 
vectors (13.101) . However, fixing the full family of coefficients is a hard task. It was performed for small TV in 
the case of the XXX chain on the segment [5], by using the a priori knowledge of the eigenvalues. It appears 
that the Bethe vector factorizes in terms of a modified creation operator. This modified creation operator 
and the associated family of other modified operators are the key step to implement the MABA without 
the knowledge of the eigenvalues. These operators are linear combinations of the ones used to construct the 
highest weight representation. This modified family of operators reduces the transfer matrix to a modified 
diagonal form (see below) and the commutation relations are similar to the diagonal boundary case. Thus, 
a part of the calculation of the MABA is similar to the usual ABA and reduces to already known cases by 
simple limit. Indeed, for the choice of the gauge parameters, 
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(3.13) a => a dl (M) = p =► p dl (M) = 

the dynamical transfer matrix (12.211) at the point too + 2 A1 takes a modified diagonal form, i.e. it depends 
only in srf and @ operators, namely 

(3.14) t(u) = t d i(u, mo + 2 M) = a(u)&f d i(u, mo + 2 M) + d(u)S> d i{u, mo + 2 M). 

Thus it becomes natural to consider 3S d i{u, m) as a modified creation operator that we will use to construct 
the Bethe vector. In the following, to apply the ABA, we will need the commutation relations 

(3.15) 3§ d i(u,m + 2 )38 d i(v,m) = £$ d i(v,m + 2 )&8 dl {u,m), 
g/ d i(u, m + 2 )38 d i(v,m) = f(u,,v)£% d i(v,m)s/ d i(u,m) + 

(3.16) g d i(u, v, m)@,u(u, m)stf d i(v, to) + w d i(u, v, m)£8 d i{u, m)3> d i(v, to), 

@ d i(u,m + 2 )3§ d i(v,m) = h{u,v)38 d i{v,m)@ d i{u,m) 

(3.17) +k d i(u, v, m)S8 d i{u, m)3> d i(v, to) + n d i(u, v, m)A3 d i(u, m)srf d i{v, to), 

obtained from the method described in [8] and with functions defined in the appendix El See also j7j for 
details. 

An important fact is that the operators of the two families are related by linear relations. In particular 
we have 


(3.18) s/ d i(u,p) = srf dr (u,m) + c(p,m) < £ dr {u, to) + b(p,m)& d i(u,p - 2) 

(3.19) S> d i{u,p) = 3> dr (u, to) - 0(u)c(p,m)^d r (u,m) - <f>(u)b(p,m)£8 d i(u,p - 2) 


with 

(3.20) 


b(p, to) = 


Ot d r 


- q m ~P 


Otdl 


c{p, to) = 


a dr - q m+p - 2 Pdi ’ q 2 - m - p a dr -Pdi' 


q m - p P dr ~ Pdi 

7 2 —m—r, 


It allows to find the action of the operators of the modified transfer matrix (12.211) on the highest weight 
vector «32D, namely 


(3.21) ^( w » m o)|^m 0 )dr =uk (u)A(u)\Q^ o ) dr + Vm o ^ d i(u,m 0 - 2)\Q,^ o ) dr , 

(3.22) ^ d i(u,mo)\n^ o ) dr = U(j){q~ 1 u~ 1 )k~((r 1 u~ 1 )A{q~ 1 u~ 1 )\Q,^ lo )dr 

(frid^Vrno^dl (n, mo 2)|Il mo ) t ; r , 

with 


(3.23) 


Pmo = b(m 0 ,m 0 ) = 


<**■ - /3diq 2m °~ 2 ’ 


Remark 3.1. These actions have an off-diagonal contribution characteristic of the MABA already pointed 

out in 13HI El- 


4. MABA 

From the modified dynamical & d i operator, used as creation operator, we construct the string of M 
operators 

(4.1) B d i(u,m,M) = 3§di (ui, to + 2(M - 1 )).. .B8 d i{u M ,m), 

and act with it on the highest weight vector (13.21) . It yields the vector, 

(4-2) I *"(«)> = Bdi(u,mo,M)\n^ o )dr. 
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Using the commutation relations (13.1513.1613.1711 . the off-diagonal actions on the highest weight vector 
(13.2113.2211 . and the functional relations (IA.8IA.9ll to simplify the term with the additional 38 d i operator, we 
obtain from the usual algebraic Bethe ansatz arguments, see [7 for more details, the off-shell action for the 
modified transfer matrix, 

M 

(4.3) t d i(u,m 0 + 2M)\'$™ 0 {u)) = Af (u, u)\A>™ 0 (u)) + ^ F(u, uf)Ef (u u u, ; )|4'^ o ({u, u. f })) 

i =1 

+f]m 0 u~ 1 c(qu)B d i({u,u},mo — 2,M+ l)|U^ 0 )d r , 

with 

(4.4) A% J (■ u , u) = (f>(u)k + (u)k~(u)A(u)f(u, u) + </>(g~ 1 it _1 )fc + ((7~ 1 7U 1 )fc _ (g _1 TU 1 )A(g _1 TU 1 )h(it, u), 

(4.5) E^(ui,Ui) = 4>(q~ 1 u~ 1 )(j)(ui)(jc + (ui)k~(ui)A(ui)f(ui,Ui) 

-k + {q~ i uf 1 )k~{q~ 1 uf 1 )A{q- 1 uf l )h{u i ,u i ) s j, 

(4.6) T) mo = <rV7m 0 -lT7m 0 - 

and where we denote the set u = {u\, ..., ««} with cardinality M, the set Ui = {ui, u 2 ,..., Uj_i, Uj+i, ..., ««}, 
where the element rq is removed, and the set { 11 , 74 } as the union of u and Ui. 

Remark 4.1. From the requirement f] mo = 0 we recover the result of jS] which implies constraints between 
left and right boundary parameters. See m for more details about the limit cases. 

The new feature, characteristic of the MABA, is the presence of the string of M + 1 creation operators in 
the off-shell equation (14.31) . This term can be handled according to the following proposition, 


Proposition 4.1. For M = N, the length of the spin chain, the string of N + 1 modified creation operators 
&di acting on the highest weight vector (3.2 1) has an off-shell structure given by, 

fj mo u~ 1 c{qu)B d i({u,u},m 0 - 2, N + l)|U^ 0 ) dr = A^(u,u)|^ 0 (u)) 


N 


(4.7) 

+ y ^F(u,Ui)E^ (yti,Ui) |T 

i =1 

with 


(4.8) 

A f («,«) = Xc('u)c(q~ 1 u~ 1 )A(u)A(q- 1 u~ 1 )G{u,u), 

(4.9) 

c ( u i)c(q~ 1 uf 1 ) . . _i _! 

E (ui,ui)= x u 2 \ A K ) A [q u { )G(Ui 

y mui) 

and 


(4.10) 

~ ~(KT . KT . N+l . ih -N-l\ 

X=-KKTT — -(- + ^q + +7T<7 

\KT KT 1 


N ( 
m 0 V 


The function A (u) is given by \3.6}) and the constant factor r) mo by (O- The auxiliary functions c(u), b{u), 
F(u,v) and G{u,v) are collected in appendix\A J see equations if A.l\ ) and (A. 21) . 


The proof of this proposition will be postponed to the next section. 


Remark 4.2. This new term gives the inhomogeneous part of the modified T-Q Baxter equation introduced 
in the ODBA approach |10j and recovered from the So V approach in [20] . 
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Remark 4.3. The expression of the off-shell action of the modified creation operators 3§di conjectured in [J] 
have additional terms. These terms come from the shift of mo by —2 of the l.h.s of the equation ra here by 
comparison of the l.h.s of the equation (5.13) in [7]. This difference is due to the way we apply the MABA: 
in \ 7}, we put the transfer matrix into a triangular modified form and here we put the transfer matrix into a 
diagonal modified form. In both cases the off-shell Bethe vector are the same, up to an overall constant. 


The proposition 14.11 allows us to write the main result of this paper, i.e., the expression of the off-shell 
action of the transfer matrix with two general boundaries on the Bethe vector (14.21) . 


( 4 . 11 ) 

with 


N 


t(u)\*«(u))=A«(u,u)\*"(u)) + y F{n, Ul )E N 


( 4 . 12 ) 


A N {u, u) = A% (u, u) + A(?(u, u), 


E N { Ui 


) = Eff (ui,uf) + Ef: (ui, uf) 


where the modified diagonal part is given by (14.4I4.5[) and the off-diagonal contribution by (14.814.91) . For 
E n ( ui, Ui) =0, which are the Bethe’s equations of the model, we obtain the on-shell action of the transfer 
matrix which describes the eigenproblem of the transfer matrix (12.181) . 


5 . PROOF OF THE OFF-SHELL ACTION ( 14 . 71 ) 

In this section we prove the off-shell equation which is the core equation of the MABA. We will use 

the left SoV basis of the dynamical 38di operator constructed in II6| . Let us consider the left-vector given 
by the product of the contravariant vectors (12.101) . 

( 5 - 1 ) di(^m | = ®n =1 di{Y(v n ,m + n) \ 

It is shown in [16] that it is a pseudo-eigenvector of the dynamical 3§di operatoiQ , namely, 

(5-2) di(n%,\&di(u,m) = dt{n%_ 2 \ri%A b (u) 


with 

( 5 . 3 ) 

and 

( 5 . 4 ) 


N 



.■f*T\ 7m 

fir ) 7™+n7™+n+i 


N 

Aft(zt) = uc(u) b(qu/v„)b(uv n ). 

n= 1 


It follows that the vectors 


( 5 . 5 ) dl(&m(h) \~ dli^m^dl 1 ( V 1 ) *71 + 2 ) . . . ^ di N (v N ,771 + 2 ) 

with h = {hi ,..., hw} and hi £ {0,1}, i = 1,..., N, are pseudo-eigenvectors of the dynamical £$di operator, 


(5.6) 


di(&m(h)\&di{u,m) 


N 


7m A b(u)Y[f( 




dl (fit m—2 ifi) | 


and its pseudo-spectrum is simple for generic inhomogeneity parameters, see m for details. Thus the SoV 
basis (15.51) describe the full Hilbert space TL. 

Then, using the left-SoV-basis (ESI), we can project the equation (1+71) on the 2 n one-dimensional subsets 
of the Hilbert space TL. It reduces the equation (14.71) to functional relations which we can prove by usual 


^The term pseudo-eigenvector rather than eigenvector is due to the shift of m. 
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analytical methods. In fact, let us first consider the projection of the string of dynamical operators 
From the properties of the SoV-basis we find, 

(5.7) 


(^mo+ 2 (iv-i) (h)\B dl ({u, u}, Too - 2, N + 1) |fi mo )d r = 

N+l N 

A b (u)A b (u)W(h) Vmo+2(N-i)(Y[f( V 7 1 ^ U ) hi f( V 7 1 ^) hi ) dl{^ 0 - 4 \^ 0 )dr, 


1=1 
N 


(5-8) dl{^mo-\-2(N— 1) (^) N) |^ mo )dr — 


N N 


i=1 


i =1 


where we have used the formula, 

(5.9) 

with 


dl (* m (h)\n"U r = W{h) dl msi") dr 


N 


(5.10) w(h) = n (^ rl * _ K rl ) A K rl )) 4 

i =1 

which follows from the recursion relatior@, 

(5.11) dl(^m(h)\£l mo )dr = (v N k ( V N )A(v N dl (^m-2 (W) | ^m 0 ) dr 


and where we have 


N 


N 


(5.12) di{^m\^m 0 )dr = di(Y(v n ,m + n)\X(v n , m 0 + ri)) dr = 


q(a dr q- m °~ n - p d iq m+n ) 


'Tm+n+1 


n=1 n =1 

By means of these projections, we reduce the matrix equation (14.71) to the functional equations, 

?N( 


(5.13) -xu 1 c(qu) = 

where we use 

(5.14) 


Ag(u,U,) , A,/ v E gi.Ui,Uj) 


YF(u,m )- 


Mu)nunv-\u) h > u ' ■ ■ / A h ( Wi )n;=i/(^ 1 -^) hj 


. iv «a<^ 0 - 4 l^ 0 )«ir 
VmoVmo-2 = “X- 


dl (^m 0 —2 l^m 0 ) dr 

Using the explicit form (14.814. 91) . we reduce the proof to the checking the 2 N functional relations, 
( 515 ) ^ 1 _ P{uj,v,h) _ P(u,v,h ) _ 1 


jr[ b(quui)b(u/ui) b(quiu n )b(ui/u n ) J]^ =1 b{quu n )b{u/u n ) 


with 

(5.16) 


N 


P{u , v, h) = ( b(uvi)b(qu/vi )) * ( b(quvi)b(u/vi )) 


i-/ii 


i=l 


It can be proved by considering left and right side of the equation (15.151) as two meromorphic functions of 
u and showing that they have the same poles and residues (at u = Ui and u = on the full complex 

plane and the same values at u = oo and u = 0. One uses the property P{u, v, h) = P(q~ 1 u~ 1 ,v, h). 


^This recursion relation is proved in m- 
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6. Connection with the ODBA and the SoV methods 


The ODBA and the SoV methods use the analytic properties of the transfer matrix. It is a Laurent series 
in u of degree 2 N + 4 and satisfy the parity condition and the crossing symmetry, see e. g. im 

(6.1) t(—u) = t(u), t(q~ 1 u~ 1 ) = t(u). 


It follows from these two properties that the transfer matrix can be rewritten as a polynomial of 

(6.2) U(u) = 

of degree N + 2 with coefficients U acting on H. 


qu 2 + q 1 u 2 
' + q~ 1 


N +2 

(6.3) t(u) = tiU(u)\ 

i =0 

We note that this variable U (u) is used to extract the q-Onsager algebra from the reflection equation [2] [3]. 
To fix any polynomial of degree N + 2 one has to find N + 3 values at distinct points. From the explicit 
form of the transfer matrix (12.41) in terms of the R-matrix (12.11) . the K-matrices (12.31) and (12.21) . it is a simple 
exercise to find that 

N 

(6.4) f(l) = (q + q~ 1 )(i ,+ + v~)(e+ + e~) Hqu^qy- 1 ), 

2=1 


(6.5) 


N 

t{i) = {q + g _ 1 )(p“ - u+ )( e ~ ~ e + ) II KiWiMkvy 1 ), 

j =i 


( 6 . 6 ) 


4+2./V 2 +N 

lim t(u) = y - -\\ 2 n (^ 2 ^ 2 + « 2 t 2 ) + 0(u 2+2N ) 

u—> oo (q — q 


which fixes 3 points. To get the remaining N points the properties of the R-matrix and the K-matrices, see 
uni for the proof, allows to relate the transfer matrix with the Sklyanin determinant from the relation 

Det q (K + (q- 1 v i )K(q- 1 v i )) 


(6.7) 

where 


t(q Vi)t(v.i) = 


b{qv^)b(qv i 2 ) 


(6.8) Det q (K + (u)K(u)) = b{u 2 )b(q 4 u 2 )k + (qu)k ( qu)k + (q l u 1 )fc (q 1 zi 1 )A(qu)A(q l u 2 ). 

This relation is the core of the SoV and the ODBA methods. Using the polynomial properties of the transfer 
matrix, the t(q~ 1 Vi) can be related to the t(Vi) by Lagrange interpolation (see below the SoV characterization 
of the spectrum). It leads to N quadratic relations on the t(vi) and gives at most 2 N solutions from the 
Bezout theorem. The ODBA proposes an analytical ansatz for the Baxter T-Q relation that satisfy the 
conditions m , m, (Esi). esi and m evaluated on some unknown eigenstate of the transfer matrix. 


6.1. Connection with the ODBA. The key step in the ODBA was to propose a modified (or inhomoge¬ 
neous) Baxter T-Q equation of the form 

(6.9) A N (u, u)Q N [u , u) = ^{u)Q n ( q~ 1 u, u) + ip(q^ 1 u~ 1 )Q N (qu, u) + xc{u)c(q~ 1 u~ 1 )A(u)A(q^ 1 u^ 1 ) 
with 

(6.10) Q n (u) = G(rt,zi) _1 , ip(u) = cj)(u)k + (u)k~ (u)A(u) 

and the Bethe equation (14.121) ensuring the analyticity of the eigenvalues. The most important point is 
the addition of an inhomogenous term in the ansatz such that the conditions (16.11) . (16.41) . (16.51) . (16.61) and 
(16.71) evaluated on some unknown eigenstate of the transfer matrix are satisfied. In the MABA context, the 
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relation (16,91) appears naturally as a consequence of the algebraic properties of the double-row monodromy 
operators, see (14.121) . 

In ^10. another ansatz depending on the parity of the length of the chain N was proposed, actually it 
appears that there is an infinite family of modified Baxter T-Q equations which parametrize the spectrum 
of the transfer matrix. Here, from the MABA we recover the minimal solution with N Bethe roots (EH) 
that is equivalent of (14.121) . This minimal modified Baxter T-Q equation and the associated Bethe equations 
were shown from the SoV to give the complete description of the spectrum for the inhomogeneous case [20] 
and numerical evidence show that it should remains the case in the homogeneous limit. 

Moreover, from the modified Baxter T-Q equation (1 6.9 [) and projection of the Bethe vector with the left 
SoV basis similar to the one we use here, the on-shell Bethe vector has been obtained in [33,134]. The MABA 
gives the algebraic construction of the off-slicll Bethe vector and allows, by limit or by imposing constraints 
between left and right boundary parameters, to recover the previous results in the literature, see |7] for 
details. 


6.2. Connection with the SoV characterization of the spectrum. The SoV representation of the 
spectral problem of the transfer matrix is given [IB] (up to some restrictions for the values of the boundary 
parameters and of the inhomogeneous parameters, see in m) by 


( 6 . 11 ) 

with 

( 6 . 12 ) 

and 

(6.13) 


AW ( U ) = J2 AN ( l ’j)9j(u) + f(u) 
1=1 


TT u ^ - U 

U ^) 2 -^ k }} k ^U{v 3 )~U{v k ) 


f(u) = A N (1) 


U(u ) + 1 

U{ 1 ) + 1 


n 


U(u) - U(vk) 
1 - U(vk) 


A AT CD U ( U ) ~ 1 TT ~ U ( U ) 

[ U(i)- lW 1 + U(v k ) 


(q + q~ 1 ) N+2 

( q-Q- 1 )™ 


(k 2 t 2 + k 2 t 2 )(U(u ) 2 


N 

fc=l 


provided that the A N (vj) satisfy the quadratic relations induced by the relation (16.71) and the interpolation 
formula (16.111) at u = g -1 ^. It is proven in j20j that this representation is complete and that it is equivalent 
to the modified T-Q Baxter equation (16.91) . 


6.3. Connection with the SoV states. To give the SoV states, let us recall the right-SoV-basis used to 
construct the eigenstates of the transfer matrix. Following m we introduce the right pseudo-vacuum 

(6-14) \^Z,)di = ®n=i\ Y ( v n,m + n)) d i 


that is a pseudo-eigenstate of the 9§di operator 

(6-15) &di(.u,m)\n*) d i = V^A b (u)\n^ +2 ) d i 

with 

(6.16) f = -QT 2 (q~ N+m Pdi - i^) (q~ N+m Pdi - i^) 

V flT/ \ flT / 7^ +1 

and 


N 

A b (u) = u c(u) b(quv n )b(u/v n ). 

n= 1 


(6.17) 
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= — — (X(u,m + 2)\K(u)\Y(u 1 ,m)) 
*7 m 


The other states are constructed from the action of the operator 
( 6 . 18 ) 

and are given by 

( 6 . 19 ) \^m(h))di = II {&di(vi,m)^ |n^)dz. 


N 


i=i 


The action of the &di(u, rri) operator on these states is 


( 6 . 20 ) 


N l _ h 

3S d i(u, m)\^f m (h)) d i = fj™A b (u) (f(vi, u)) \^ m+2 {h)) M . 


i= 1 


The SoV states are given in m , up to an overall normalization, by 


( 6 . 21 ) 


N 


A N (V~ 1 ) y.h 


^ ^(/j) h[ { Vj k+( V . ^(y. i) 


4-rprrr) 3 1 * m+ 2 Ch)U, 

\(D\ V- ]/ 


where the A N (Vj 1 ) are fixed from the quadratic relation (16.71) and the interpolation formula (16.111) . We note 
that in m the gauge parameter j3 is arbitrary and we restrict here to the case /3 = Pdi to make the link 
with the Bethe vector (14.21) with M = N. The /i„(/i) is the diagonal element of the scalar product between 
left and right SoV basis 

N 

( 6 - 22 ) di(* m (h)\* m +2(k))di = V*(h) 


explicitly given by 


(6.23) 


where 


.N 


Mr, 


(j)=n 


N 


dl 


i=l ^m+l+i 


( Vi)b{q 


-2hi+l v 2 


1) 


N 


hhq 

j<i 


-2hj+l 


Vj/Vi)b(q 


- h i- h <+\/v.)b(q-* h i+\ Vj )b{q h i- h * 


J i v 3) 


(6.24) 


rydl rydl 

-N N Im+N /m+l+AT 
Vm=Vm- 


% 


dl 


To relate this result with the on-shell Bethe vector m with M = N constructed from the MABA we 
project both vectors with the left SoV-basis. For the Bethe vector we use the formula (15.81) with u = Ui 
and for the SoV vector we use the scalar product between left and right SoV-basis (|6.22ll . Thus we obtain 
in both cases, up to a overall constant, the same function, for each of the 2 N components, that show the 
equivalence of the two characterizations for the eigenstates of the transfer matrix. 


7 . CONCLUSION 

In this paper we have proved the expression of the off-shell action of the modified creation operator on 
the Bethe vector, see nroDQsition l4.il Thus we give the proof of the off-shell action of the transfer matrix of 
the XXZ spin chain on the segment conjectured in (T) and before in some limit cases 03®. The proof uses 
the left SoV basis m and exhibits similarities with the proof of the on-shell Bethe vector from the ODBA 
method [34] . 

This shows the deep relation between the MABA, the ODBA and the SoV methods. Indeed, the MABA 
explains the algebraic origin of the inhomogenous term in the modified T-Q Baxter equation introduced in 
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[ID] from the ODBA, as a modification of the analytical Bethe ansatz. Moreover it provides a correspondence 
between the Bethe vector and the eigenvectors in the SoV framework, completing the link already done for 
the eigenvalues in !2DJ. Let us add that it implies the existence of a modified coordinate Bethe ansatz whose 
wave function can be obtained by projection of the on-shell Bethe vector on the appropriate basis, extending 
the constrained boundary cases obtained in Ean2 to generic boundaries. Thus the MABA incorporates 
the different BA (ODBA, Analytical BA, ABA, Coordinate BA) and is equivalent in the inhomogeous case 
to the SoV. 

The SoV method yields a proof of the completeness of the spectrum, at least in the inhomogenous case, 
which can be extended to the MABA and, by complementarity, the MABA provides a regularization for the 
SoV such that the homogenous limit is well defined. We remark that such a connection was already used in 
the case of the XXX spin chain on the circle in [23], where the completeness of the Bethe ansatz solution was 
proved. Thus it seems that the complete description of the spectral problem of finite dimensional quantum 
integrable models constructed from the Yang-Baxter and reflection equations can be handled by both the 
Bethe ansatz and SoV methods. 

It could be also of interest to relate the MABA for the XXZ spin chain on the segment with other known 
methods such as the q-Onsager approach [2] or the non-polynomial solution from the homogeneous Baxter 
T-Q relation [2T]. 

The knowledge of the Bethe vector is the first step to consider the correlation functions of the XXZ spin 
chain on the segment in the ABA framework. An important question will be to understand the structure 
of the scalar product of the Bethe vector and see if some determinant representation can be obtained by 
analogy with the diagonal case [19]. Let us mention the recent algebraic functional approach [18, used in 
the diagonal case that should allow to find integral representations of the scalar product from our off-shell 
results. 

Other interesting questions will be to apply the MABA to other models such as periodic spin chains 
without U( 1) symmetry where the ODBA [9] [33] and the SoV _31, 27] methods were applied as well as 
developing a Nested MABA for spin chains without 17(1) related to higher rank algebras where only the 
spectrum was obtained in few cases from the ODBA, see e. g. El- 
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Appendix A. Functions 
We use the following functions throughout the text, 

-l 


(AT) 

(A.2) 

(A.3) 

(A.4) 


b(u) = 

<t >( u ) = 


u — u 


-—, k ( u ) = v-U + v + u 1 , k + {u) = e+u + e-U 1 , c[u ) = u 2 — u 

q~q 

b(q 2 u 2 ) 


f(u,v) = 


b(qu 2 ) 
b(qv / u)b{uv 


G(u, v ) = 


b(u/v)b(quv) 

(friq-'v- 1 ) 


F(u,v) = (v/u)G(u,v) 


b(q 2 u 2 ) 

4>{v) 


b(v / u)b(quv) 


, g(u,v) = 


*<»,«) = 

b{quv)b(u/v) 


b(u/v) 

/>(■») 
b(v/u ) ’ 


, w(u,v) = - 
n(u, v) = 


1 


b(quv ) ’ 
(f>(u)cj)(q~ 1 v^ 1 ) 
b(quv) 
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(A.5) 

(A. 6 ) 


7 (u,m) = aq rn u-/3q m u 1 , j(l,m)='y m , 

r i(u/v,m + 1 ) j(uv, to) 

g{u,v,m) = - g(u,v), w(u,v,m) = - w(u,v), 


7m+l 


7m+l 


,. ,, , j(v/u,m+ 1 ), . . , . 7 ( 1 /(«v), in + 2 ) 

(A.7) k(u,v,m) = - -k(u,v), n(u,v,m) =- -n(u,v), 


7m+l 

The following functional relations hold 

M 


7m+l 


(A.8) /(u, u) + ^ s(u, Ui, m - 2)/(«,, Uj) - <j>(v,i)w(u, Ui,m — 2)h(ui , u,) = 


i=l 


7m-2M-l 

7m-l 


and 


M 


(A.9) ft(u, u) + <£(u) 1 ^ <f>(ui)k(u, Ui,m — 2)h(m, m) - n(u, m, m - 2 )/(«*, u,) = 


7m-2M-l 

^m— 1 


To prove these functional relations we consider them as a function of u. One can easily check that the left 
hand side of these equations are holomorphic functions on the full complex plane, thus they are constant. 
To evaluate these constants one sends u to infinity and use the relation 


(A.10) 


M 

i —1 


q 2 M - q~ 2M 

n — n 1 


This relation can be proved in a similar way. The left hand side is a symmetric function of the set iZ, so one 
can consider it as a single variable function of any Uk and show that for generic parameter Uk this function is 
holomorphic on the full complex plane, thus it is constant. Then sending Uk to infinity one find a recursion 
relation on M that fix the constant. 
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